Topologically protected surface Majorana arcs and bulk Weyl fermions in 

ferromagnetic superconductors 
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A number of ferromagnetic superconductors have been recently discovered which are believed to 
be in the so-called "equal spin pairing" (ESP) state. In the ESP state the Cooper pairs condense 
forming order parameters A^^,A|| which are decoupled in the spin-sector. We show that these 
three-dimensional systems should generically support topologically protected surface Majorana arcs 
and bulk Weyl fermions as gapless excitations. Similar protected low-energy exotic quasiparticles 
should also appear in the recently discovered non-centrosymmteric superconductors in the presence 
of a Zeeman field. The protected surface arcs can be probed by angle-resolved photoemission 
(ARPES) as well as fourier transform scanning tunneling spectroscopty (FT-STS) experiments. 
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Introduction: Ferromagnetic superconductors (FS) 
display a remarkable coexistence of the mutually ex- 
clusive order parameters of ferromagnetism and super- 
conductivity in the same band of electrons [H, In 
the recently discovered Uranium-based Ising ferromag- 
netic superconductors UGe2 @, URhGe [l], and UCoGe 
0, it is believed that superconductivity appears in the 
spin-triplet p-wave equal-spin-pairing channel Q . In this 
channel the superconducting order parameters A-|-|- and 
A4,4, are decoupled in the spin space. In FS typically the 
ferromagnetic transition temperature T* far exceeds the 
superconducting transition temperature Tc, and the com- 
plex interplay of the spin-triplet superconductivity and 
ferromagnetism leads to a fascinating phase diagram as 
a function of temperature, pressure, and magnetic field. 
In this paper we show that in 3D the superconductivity 
in these materials generically supports protected chiral 
Dirac or Weyl nodes in the bulk and topologically robust 
open Majorana fermion arcs on the surface. Similar bulk 
and surface gapless quasiparticles, but without the Ma- 
jorana character, have been recently predicted in some 
non-superconducting systems . They have also been 
predicted to occur in superfluid He3-A and the re- 

cently discovered inversion-asymmetric superconductors 
known as non-centrosymmetric superconductors (NCS) 
12 -[l4 1 in the presence of time reversal (TR) invariance 
15|. As we show below, the FS systems are another can- 



didate providing one of the simplest platforms for the 
realization and detection of 3D gapless topological su- 
perconductivity (TS) with protected surface Majorana 
quasiparticles. Such modes on the surface should be 
accessible to surface sensitive probes Fourier transform 
scanning tunneling spectroscopy (FT-STS) [iMi- 

The Weyl nodes in the FS are characterized by an en- 
ergy dispersion linearly proportional to momentum and 
come with a specific handedness or chirality. There are 
an even number of such nodes with each pair consisting 
of nodes with opposite chiralities. Such chiral Dirac, or 
Weyl, nodes are topologically protected in 3D by an in- 



variant [19] which takes the values ±1, and hence can 
only be removed when a pair of such nodes with oppo- 
site signs of the invariant collide in the momentum space. 
From the bulk-boundary correspondence such topologi- 
cally protected gapless bulk spectrum leads to open Ma- 
jorana fermion arcs on the surface which should be de- 
tectable in STM tunneling experiments. 

Recently, there has been a lot of excitement about 
two- and one-dimensional TS systems with broken TR 
symmetry which can be used to realize interesting non- 
Abelian quasiparticles such as Majorana fermions. Zero- 
energy Majorana fermions, defined by hermitian opera- 
tors 7g — 7o, can in principle be used to realize fault- 
tolerant architecture for quantum computation [g^l . The 
canonical example of these systems is the 2D spin-less 
Px + ipy superconductor [lH which is gapped in the bulk 
and has a single zero-energy Majorana fermion mode lo- 
calized at order parameter defects such as vortices. Addi- 
tionally, this system has gapless chiral Majorana fermion 
modes (defined by operators 7]^ = 7-fe) on the sample 
boundary. The robustness of such non-trivial topolog- 
ical properties of a 2D spin-less Px + ipy superconduc- 
tor can be understood in terms of the existence of a Z2 
topological invariant which can be viewed as the parity 
of the Chern number of the corresponding 2D particle- 
hole symmetric Bogoliubov-de Gennes (BdG) Hamilto- 
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23j . Practical realizations of these exotic prop- 
erties have recently been proposed in strong topological 
insulators (TI) which have a strong spin-orbit coupling 



24| and, remarkably, even in ordinary semiconductors 



with a strong spin-orbit coupling and a suitably directed 
Zeeman splitting [25l - E7j |. 

In this paper we extend the concepts of TR-breaking 
TS states to 3D solid state materials and investigate pos- 
sible experimentally accessible condensed matter systems 
where such physics can be realized. We find that FS 
systems are ideal candidates in 3D where very similar 
physics is realized even in the absence of a spin- orbit cou- 
pling and an external Zeeman field. Interestingly, we find 
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that in FS the analogous states are no-longer gapped in 
the bulk, but have gapless points in the momentum space 
with a Weyl spectrum of the BdG quasiparticles. Even 
more interestingly, the existence of the bulk Weyl nodes 
directly corresponds to the existence of open gapless Ma- 
jorana fermion arcs on some suitable surfaces of the 3D 
system. The open surface Majorana fermion arcs offer 
the tantalizing possibility of detecting gapless Majorana 
excitations using the available surface sensitive probes 
such as ARPES and STM tunneling experiments. 

Topological phases in ferromagnetic superconductors: 
Ferromagnetic superconductors, unlike conventional s- 
wave superconductors, are characterized by a spin-triplet 
pairing potential which separates in the spin-sector (i.e., 
A-|-i- and A^j, are the relevant order parameters). More- 
over, because of the existence of a strong internal TR- 
breaking exchange field far greater than the usual Pauli 
limit, the superconducting order is thought to be [i^ of 
the so-called non- unitary p-wave type Q . Based on these 
constraints, we will take below the relevant representa- 
tive order parameter as A^^i{k) — Saa'^ir^^^^^^^, with 
^ , which is the order parameter of the A2 phase of 
He3 . The particular orbital form of the order param- 
eter, however, will only be used at the end for numerical 
calculations. Most of the general results below (includ- 
ing the important predictions of the bulk Weyl modes 
and surface Majorana arcs) follow more generally simply 
from the equal-spin-pairing structure of the superconduc- 
tivity of these materials. This is important because so far 
there is no consensus on the orbital symmetry of the su- 
perconducting order in FS. 

The mean-field Hamiltonian describing an ESP state 
for a p-wave superconductor is written as 

Hbcs = Y.j d'H^^^^^ ~ £F,.)/],/.fc 
+ [A^.'{k)fljl_^ + h.c] (1) 

where a =t, ^^ labels the spin index of the electron opera- 
tors and the pairing potential is odd in momentum space 
so that Acrcr' (fc) ~ —Aa'a{~k). The magnetization of the 
ferromagnetic superconductor is accounted for by the dif- 
ference in the Fermi energies (£F,t ~ ^f,])- Defining the 

Nambu spinor *(fc) = (V'|(fc), '0|(fc), V't(-'^)> ^il"'^))' 
the BdG Hamiltonian for the ferromagnetic SC is writ- 
ten as 

Hb{k,,ky, k,) = ( " - SF.av - MN{0)a,)T, 

+ [A(fc)r+ + h.c] (2) 

where M — ^^2jV(o)''^ proportional to the magnetiza- 
tion, iV(0) is the density of states at the fermi level, and 
£F.av = ('^-^■t+^-F.i) jg i^j^g average fermi energy of the two 
spin components. 



The properties of the three-dimensional phase de- 
scribed by Eq. [2] can be understood by dimensional re- 
duction in momentum space to classes of two-dimensional 
topological phases. The idea of dimensional reduction is 
to consider the Hamiltonian Hii{kx, ky, kz) for the trans- 
lationally invariant 3D bulk system as a set of 2D systems 
i.e., 

"^i^x^ky) — H()[kx, ky , kz) (3) 

where H^'^\kx,ky) are effective 2D Hamiltonians 
parametrized by kz- The topological properties of such 
2D Hamiltonians, which break the TR symmetry, are 
characterized by the Chern number associated with the 
occupied states [2^ . The Chern number of a 2D topolog- 
ical Hamiltonian counts the number of chiral edge states 
at a given kz- The edge states of the 2D system can be 
understood as surface states of the 3D system for surfaces 
which are parallel to the z-axis. For such surface states, 
the wave- vector along z, kz is a good quantum number 
and the dispersion of the surface states can be obtained 
from the dispersion of the edge states. 

Pfaffian topological invariants: The 3D BdG Hamil- 
tonian Hf) of the FS systems has a particle-hole symme- 
try A = iTyK, which anti-commutes with Hh- Here K is 
the complex conjugation operator. To make use of this 
symmetry, we will assume from here onwards that Hb is 
an even function of kz so that H^^^ is also particle-hole 
symmetric. In this case, one can use the Pfaffian topolog- 
ical invariant, sgn{P f {iTyHj^^\kj: — ky — 0))) jl^ to 
determine the parity of the Chern number, which must 
change whenever the Pfaffian changes sign. Moreover, 
since the pairing potentials vanishes at (k^ = ky = 0,kz), 
the Pfaffian topological invariant is found to be 

sgii{P f{iTyHj^l°\kx,ky = 0))) = sgn(eKt-A:^)sgn(£F,;-fc^). 

(4) 

The above topological invariant is non-trivial only in the 
restricted range of values of kz which satisfy: 

kz,c,i = \J1m*eF,\. <\kz\< kz^c.^ = ^/2rn*e^. (5) 

For values of kz outside this range, the system has even 
Chern parity. 

Bulk Weyl fermions: The family of topological super- 
conductors described by the Hamiltonian Hj^^^ under- 
goes a quantum phase transition from the topological to 
the non-topological phase when the Chern parity changes 
at the values of kz where the conditions in Eq. [5] are 
saturated. Such topological quantum phase transitions 
are accompanied by a closing of the topological gap at 
kx = ky ^ 0- The set of energy eigenvalues of the 3D 
Hamiltonian Hf){kx, ky,kz) in Eq. [2] is a union of the en- 
ergy eigenvalues of the entire family of 2D topological 
superconductors described by the Hamiltonians H^^\ 
Therefore the Hamiltonian Hb{kx, ky,kz) must have gap- 
less points ai K = {kx = 0,ky = 0,kz = ikz.ca), with 
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a two-fold degeneracy of eigenstates, \tz = ±l,cr2 = cr), 
where a =t, i- The dispersion of the pair of degenerate 
states { \tz = ±l,az = cr)) around the degeneracy points 
K can be obtained by expanding the Hamiltonian in 
Eq. Oas k — K + Sk to Hnear order in 5k. The resulting 
Hamiltonian takes the form of a three-dimensional Dirac 
cone, 



HbiK + Sk) = Sk ■ VAr{K)t.^ 



Sk-VAi{K)Ty 



5K. 



o{Sk'), 



(6) 



where Ajij{K) are the real and imaginary parts of A{K) 
in Eq. [21 Such Dirac cones in 3D are protected because 
continuous deformations of the bulk Hamiltonian Hi,{k), 
which can only redefine dk in Eq. [Bl cannot gap out such 
a Dirac cone. A perturbation to Eq. [2]can only get rid of 
Dirac cones from the spectrum by shifting and merging 
them in pairs [19]. This kind of protection of Dirac cones 
in D = 3 can be further represented by associating it with 
a topological invariant 19]. Bulk semi-metals in three 
dimensions with such Dirac-like point fcrmi surfaces are 
referred to as Weyl semimetals. Thus iJf, in Eq. [2] rep- 
resents a nodal superconductor which is the analog of a 
Weyl semimetal with four isolated Dirac cones. 

Surface Majorana arcs: The 2D topological supercon- 
ductor Hamiltonians H^^\ka:,ky), which are parame- 
terized by kz, have an odd Chern number in the range 
Eq. [5]and are therefore characterized by chiral Majorana 
edge states that are confined to the edge of the system. 
A surface along the x — z plane for the original 3D Hamil- 
tonian Eq. [5] is translationally invariant along the z and 
X directions and therefore has well-defined k^ and k^ mo- 
menta. Therefore the surface state with a fixed k^ and 
kx of the 3D Hamiltonian Eq. [2] is identical to the edge 
state with momentum k^ of the 2D Hamiltonian in Eq. [3] 
with kz as a parameter value. The energy of such a chiral 
Majorana mode vanishes for fc^ = at any value of kz in 
the range in Eq. [5] and therefore appears on the surface 
ARPES spectrum as a Majorana arc. 

The typical extent in kz of such a non-degenerate Ma- 
jorana arc is limited to the range given in Eq. [S] It is 
important to note that in contrast to TR breaking TS 
states in D = 2, 1, where edge Majorana fermions appear 
only in a restricted range of parameters (for example, for 
chemical potential smaller than the Zeeman coupling in 
^51427)). in 3D FS systems Majorana arcs appear 
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for any value of the chemical potential (albeit in a re- 
stricted range in momentum space, Eq. [SJ. Below, as 
mentioned earlier, we take the order parameter of the FS 
systems to be of the non-unitary ESP type with an or- 
bital structure given by Aa-a-'{k) ~ Saa'Aa- with 
A-f 7^ A^. In this case, the pairing potential A^^'{k) is 
associated with a Chern number 2, therefore one would 
have a pair of chiral surface modes propagating along the 
surface in the range \kz\ < fcz,4,- This could in general ap- 
pear as a pair of fermi arcs of Bogolibov quasiparticles in 



the ARPES or STM spectrum. Since the BdG Hamilto- 
nian now decouples into a spin-'|" and spin- J, sector, the 
2D Hamiltonian H^^J{kx,ky) for each spin a can be 
thought of as independent odd Chern number topolog- 
ical superconductors which has a Majorana chiral edge 
mode in the ranges 



for spin-up electrons and 



(7) 



(8) 



for spin-down electrons. Therefore we find that the sur- 
face arcs exist over a much larger range in kz i.e between 
-kz^ct s-iid kz,c,t- 

These chiral edge modes exist in the spin-cr sector for 
kz satisfying \kz^a-\ < kz,c,a and have a dispersion of the 
form ek^{kx) = v{kz)kx- Therefore the dispersion of the 
corresponding surface mode is given by 

A,{kz) 



e{kx,kz,a) = v{kz,a)kx 



^2m*(, 



£F.cr 



(9) 



for \kz\ < kz^cM and 

\e{kx,kz,<j)\> eg (10) 

where eg is a finite positive gap for \kz\ > kz^c.a- The 
dispersion of the surface modes given in Eq. [9] has the 
special property that the energy vanishes on a pair of 
lines kx = Q which terminates at kz = ±kz^c,a- These 
lines are referred to as Majorana arcs. 

As mentioned before the Majorana character of the 
surface arcs is only protected in the restricted range in 
Eq. [S] In the rest of the kz range, the dispersion arcs 
are not strictly speaking Majorana in character. In par- 
ticular an in-plane Zeeman splitting along x, Vz^x, can 
lead to a mixing of the a =t, 4- states so that the surface 
Majorana arcs now split into a pair of fermi arcs, indexed 
by s = ±1, with dispersion 

(I. I. \ v{kz,^) + v(kzA) ,^ 

eyKx, Kz7 ^ - Kx 




v{kzA)~v{kzA) 



1.2 



(11) 



for \kz\ < kz,c,i- However, one must note that while 
these modes are not Majorana in the strict sense, they 
still have equal weights of electron and hole states. 

The Majorana arcs are obtained by solving 
£{kx,kz, s) = 0. This leads to the equation for the 
Majorana arc 



±- 



Vz,x 



± 



Vz.x^/ep - kl 



(12) 



' ^v{kzA)v{kz,^) ^TA^^ 

where we have assumed the magnetization to be small 
compared to the total density (i.e. — ej,] <C SF,av) and 
the expression Aa{kz) — A^. The resulting Majorana 
arcs are ellipses in the {kx, kz) plane between \kz \ < kz.d 
as plotted in Fig. 1. 
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FIG. 1: yz,a;-dependent surface fermi contour plotted ac- 
cording to Eq. [l2]for A4, = 0.6At,At = Aq. The width of 
the contour in the kx direction is proportional to Vz.x- As 
Vz,x — > 0, the spin-t and spin- 4, sectors decouple and the 
fermi contour evolves into a pair of surface Majorana arcs. 



Quasiparticle interference of Majorana arcs in STM: 
The surface Majorana arcs shown in Fig. 1 should in 
principle be directly visible as arcs in momentum space 
at the Fermi energy in the ARPES spectrum. However, 
since these arcs are separated from the bulk states by a 
relatively small superconducting gap (~ IK 0.1 meV), 
it is not clear if the typical energy resolution achieved 
in ARPES is sufficient to resolve the Majorana arcs. On 
the other hand, low temperature and low noise STM mea- 
surements often have sub-100 mK energy resolution al- 
lowing one access to the energy scale of the Majorana 
arcs. While conventional STM, which measures the local 
quasiparticle density of states at the surface, would be 
able to detect the presence of Majorana arcs as gapless 
modes at the surface, it does not provide any information 
about the structure (finite extent in momentum, curva- 
ture, etc.) of the Majorana arcs in momentum space. 
Such information can be obtained by FT-STS 16- 18]. 
FT-STS relies on the fact that impurity-scattering at 
the surface leads to a spatially varying local quasiparticle 
density of states, n{r) at the surface, which can be deter- 
mined from the spatial variation of the tunneling current 
/t(r) oc n(r) at the surface. The resulting current map 
obtained from STM can be Fourier transformed to obtain 
It{q) — J drlt{r)e^''''^ oc / drn{r)e'''^''^ . The disorder av- 
eraged square of the Fourier transform (|/(q)p) can be 
shown to be related to the joint density of states p(q) 
(i.e. p{q) cx (|/(q)P)) defined by 



p{q) = j dkS{e{k))S{e{k + q)), 



(13) 



where e(fc) is the surface mode dispersion in Eq. 11. 
The application of the FT-STS method [SM] outlined 
above for the weakly disordered surface of a ferromag- 
netic semiconductor is expected to lead to the character- 
istic structure that is plotted in Fig. [2] 



-0.25 



FIG. 2: (color online) The joint density of states p{q) (in 
arbitrary units) at the fermi surface, which, as explained in 
the text, would be measured by FT-STS on the surfaces of 
FS systems. We have used Eq. 1131 and the parameter values 
Vz,x = O.lA-f and Aj, = 0.6A; (i.e. same as Fig.[l|. 



Summary and Conclusion: We propose the recently 
discovered ferromagnetic superconductors [iHl] as exper- 
imentally accessible 3D systems where the physics of TR- 
breaking TS states is realized. These systems support 
robust bulk gapless quasiparticles in the form of chiral 
Weyl fermions and protected open Majorana fermion arcs 
on suitably oriented surfaces of the 3D system. In con- 
trast to the TR-breaking TS states in D = 2,1 [13^ 



where Majorana modes appear only in a restricted range 
of the chemical potential smaller than an applied Zee- 
man field, remarkably, in the 3D systems considered here 
edge Majorana modes appear for all values of the chem- 
ical potential. The Weyl nodes in the bulk are topolog- 
ically protected because they separate momentum space 
regions where the BdG Hamiltonians have different val- 
ues of a relevant topological invariant, the analog of the 
2D Chern number. Therefore the Weyl nodes arise from 
topologically unavoidable closing of the quasiparticle gap 
at isolated points in the momentum space. The existence 
of the bulk Weyl nodes directly corresponds to the exis- 
tence of open gapless Majorana fermion arcs (Fig. 1) on 
suitable surfaces of the 3D system. The surface Majorana 
fermion arcs offer the tantalizing possibility of detecting 
gapless Majorana excitations using the available surface 
sensitive probes such as ARPES and spectroscopic STM 
experiments (Fig. 2). 

The FS systems are not the only materials which can 
support 3D gapless TS states with broken T R sy mmetry. 
The newly discovered 3D NCS materials 12l4l4j can also 
support such states in the presence of a sufficiently large 
Zeeman splitting. In NCS (e.g., LisPd^Ptg-^jB, Y2C3), 
because of an intrinsic spin-orbit coupling, the supercon- 
ducting order parameter has an admixture of s-wave and 
p-wave components. In the presence of a Zeeman split- 
ting larger than the s-wave component of the order pa- 
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rameter, the NCS enters into a TR-breaking TS state 
with bulk Weyl fermions and surface Majorana arcs. The 
Zeeman splitting, however, should not be accompanied 
by a large orbital depairing field which may destroy the 
superconductivity itself. This can be ensured by choosing 
materials with a large enough g-factor so that a relatively 
small magnetic field can still create a Zeeman splitting 
larger than the s-wave part of the order parameter. 

J.S. thanks the Harvard Quantum Optics Center for 
support. S.T. thanks DARPA and NSF for support. 
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